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Abstract 

We study the time dependence of the ionization probability of Rydberg atoms 



in 
^t- 

\o 
o 

| driven by a microwave field. The quantum survival probability follows the classi- 

cal one up to the Heisenberg time and then decays inversely proportional to time, 
due to tunneling and localization effects. We provide parameter values which 
should allow one to observe such decay in laboratory experiments. Relations to 
the 1/f noise are also discussed. 

The pioneering experiment of Bayfield and Koch performed in 1974 Q attracted 
a great interest to ionization of highly excited hydrogen and Rydberg atoms in a mi- 
crowave field. The main reason of this interest is due to the fact that such ionization re- 
quires the absorption of a large number of photons (about 20 — 70) and can be explained 
only as a result of the appearence of dynamical chaos and diffusive energy excitation 
in the corresponding classical system. Such classical diffusive ionization requires many 
microwave periods and quantum interference effects can suppress this diffusion leading 
to quantum localization of chaos [[J. 

More recent experiments || with alkali Rydberg atoms showed an algebraic time 
dependence of the survival probability: P(t) oc t~ a , with a ~ 0.5. This result cannot be 
explained within the picture of diffusive ionization in the domain of classical chaos. The 
origin of the slow algebraic decay was attributed to the underlying structure of classical 
mixed phase space composed by integrable islands surrounded by chaotic components: 
Chaotic trajectories can be trapped in the vicinity of the hierarchy of regular islands 
M and this slows down the ionization process. 

Recent studies of quantum Poincare recurrences for the standard map model in the 
semiclassical regime || showed that quantum P(t) follows the classical decay up to 
the Heisenberg time which is determined by inverse level spacings. For t > tjj, the 
quantum survival probability starts to decay inversely proportional to time (a = 1) and 
becomes much larger than the classical one. The power a = 1 is due to the exponentially 



broad distribution of time scales to escape from some region of phase space || . These 
exponentially broad escape times are originated by tunneling from classically integrable 
regions or by the exponential quantum localization. 

Microwave ionization of Rydberg atoms constitutes a unique possibility to experi- 
mentally observe quantum Poincare recurrences |J. Indeed the highly excited states of 
hydrogen atom in a microwave field can be approximately described by the Kepler map 
which is very similar to the standard map and therefore one would expect the same 
behavior for the time dependence of the survival probability. 

In order to investigate the probability decay for the hydrogen atom in a microwave 
field, we choose the initial state with principal quantum number no and numerically 
study the survival probability P(t) in a linearly polarized monochromatic electric field. 
The Hamiltonian reads 

v 2 1 

H=- + ezsm(ut), (1) 

2 r 

wher e and u are the strength and frequency of the microwave field, measured in atomic 
units. The quantum evolution is numerically simulated by the one-dimensional (Id) 
model of a hydrogen atom and by the 3d model for atoms initially prepared in states 
extended along the field direction and with magnetic quantum number m = 0. In the 
Id model the motion is assumed to take place along the field direction (z-axis, with 
z > 0). 

In order to compare classical and quantum dynamics it is convenient to use the 
scaled field strength e = en$ and frequency Uq = urn^, which completely determine 
the classical dynamics. The classical limit corresponds to H e s = h/no — > 0, at constant 
6q, Uq. In classical mechanics diffusive ionization takes place for fields above the chaos 
border: eo > e c ~ 1/(4:9lj^ 3 ) @. Quantum interference effects can suppress this 
diffusion leading to quantum localization of chaos 0. Such dynamical localization 
leads to a quantum probability distribution fx exponentially localized in the number 
of absorbed photons N$ = (E — Eq)/uj (E electron energy, E = —1/2tIq): fx oc 
exp(— 2|A^,|/£^), with localization length 1$ = 3.3elu> 10 ^ 3 nQ 0. 

We introduce an absorption border for levels with n > n c ; like this the number of 
photons required to ionize the atom is N c = (n /2uj )(l — n 2 ) /nl). Such border occurs in 
real laboratory experiments, for example as a consequence of unavoidable stray electric 
fields experienced by the Rydberg atoms during their interaction with the microwave. 
The absorption border n c can be varied in a controlled way via a static electric field e s , 
the static field ionization border being e s n\ »s 0.13. 

In Fig. [l| we show a realistic case (no = 60, eo = 0.1, Uq = 2.6) in which, initially, 
classical and quantum probabilities decay in a very similar way and only after approx- 
imately 5 x 10 2 microwave periods the quantum survival probability starts to decay 
more slowly (P(t) oc 1/t) than the classical one which decays as l/t a , with a ~ 2.15. 
The comparison of quantum simulations for the Id hydrogen atom model and the 3d 
dynamics is shown in the inset of Fig. 0. It demonstrates that both dynamics give 
very close results, confirming that the essential physics is captured by the Id model. 
Actually, due to Coulomb degeneracy, the slow motion in the orbital momentum / acts 
as an adiabatic perturbation on the n motion and as a result the excitation in n is 
well described by the Id model 0. We put the absorption border near the initial state 
(n c = 64) in order to have p c — l</,/N c « 3.5 > 1. In this way the probability can go 
out very easily and the 1/t probability decay is observed after a short transient time of 
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Figure 1: Survival probability P(t) as a function of the interaction time t (in units of mi- 
crowave periods) for e$ = 0.1, ujq = 2.6, no = 60, n c = 64: quantum (solid line) and classical 
(dashed line, ensemble of 3 x 10 6 trajectories) Id hydrogen model. The straight lines have 
slopes 1 and 2.15, the latter coming from a fit of the classical decay for 5 x 10 2 < t < 2 x 10 4 . 
Inset: quantum survival probability for the Id model (solid line) and the 3d hydrogen atom 
(dot-dashed line). 

the order of 20 microwave periods. On the contrary, when p c < 1, strong fluctuations 
around the 1/t decay take place ||. This is analogous to the huge (log-normally dis- 
tributed) conductance fluctuations in a disorder solid with localization length smaller 
than the sample size. 

In order to confirm that the algebraic probability decay is related to localization 
effects and to the sticking of classical trajectories and of quantum probability near the 
integrable islands in the phase space, we show in Fig. |2| the time evolution of the 
survival probability distribution in the phase space of action-angle variables (n, 9) for 
the Id model. In the classical case 3 x 10 6 orbits are initially homogenously distributed 
in the angle 9 on the line n = no = 60, corresponding to the initial quantum state 
with principal quantum number uq = 60. After 50 microwave periods, the classical 
distribution of non ionized orbits shows a fractal structure which surrounds the stability 
islands (Fig. fj top left). At larger times this distribution approaches more and more 
closely the boundary critical invariant curves (Fig. Q top right). One of them confines 
the motion in the region with n > rib ~ ^o^cAo) 1 ^ 5 ~ 41 where rib determines the 
classical chaos border for given eo- Other invariant curves mark the critical boundaries 
around internal stability islands (for example at n ~ 55, corresponding to un 3 « 2). 
This phase space metamorphosis is associated with a change in the exponent of the 
power law decay of classical Poincare recurrences (the crossover occurs at t ~ 5 x 10 2 ). 
In the quantum case the value of h e s is not sufficiently small to resolve the fractal 
structure at small scales. However, the Husimi function (constructed in a way similar 
to Ref. [0]) shows similarities with the classical probability distribution at t = 50 



Figure 2: Classical density plot (top) and Husimi function (middle and bottom) in action- 
angle variables (n,9), with 30 < n < 63.5 (vertical axis) and < 8 < 2tt (horizontal axis), 
for the Id model in the case of Fig. [|. Top left: 50 < t < 60; top right: 2 x 10 3 < t < 10 4 ; 
Husimi functions are averaged in a small time interval 5t = 10 to decrease fluctuations, near 
t = 50 (middle left), t = 10 3 (middle right), t = 10 4 (bottom left) and t = 10 5 (bottom right). 
Bright regions correspond to maximal density. 




Figure 3: Quantum (solid line) and classical (dashed line) average action < n > as a function 
of the interaction time t for the case of Fig. [l]. The straight lines indicate the position of the 
stability islands for wn 3 ~ 1 (n ~ 55) and wu 3 ~ 2 (n ~ 44). 

(Fig. middle left). At longer times, the diffusion towards the boundary at n b is 
slowed down due to localization effects and penetration of the quantum probability 
inside the classical integrable islands. At t = 10 3 (Fig. middle right) the quantum 
probability is concentrated in a layer around n « 49. Due to localization effects, the 
Husimi function does not change significantly for a very long interaction time (compare 
with Fig. bottom left at t = fO 4 ). Eventually the probability starts to penetrate 
very slowly inside the main island at n « n b (see Fig. |2| bottom right at t = 10 5 ). 
Therefore tunneling and localization effects are responsible for the slow 1/t decay of 
quantum survival probability seen in Fig. |l|. The chaos border starts to influence the 
dynamics only after a very long interaction time because its distance from the starting 
line n = n is much greater than the localization length p b = i^/lN^ « 0.35 < 1, 
where N b = (n /2u )(l - nl/n 2 b ). 

Fig. m shows the average action < n > of the non-ionized electrons as a function 
of the interaction time. The diffusion towards the chaos border is slowed down in the 
quantum case and different regimes are clearly visible: tunneling inside the islands for 
oon 3 w 2 (for t < fO 2 ), localization of the wave packet (for 10 3 < t < 3 x fO 4 ) and 
tunneling inside the main island at um 3 w 1 (for t > fO 5 ). 

Notice that the probability decay P(t) is related to correlations decay via C(t) oc 
tP(t) |4[]. In the case of a = 1 this implies that correlations do not decay. The Fourier 
transform of C(t) gives the spectral density S(lu) of the effective noise produced by the 
dynamics: 

S(w) = J C(t)exp(iujt)dt oc l/u. (2) 

This shows that the spectral noise associated with the quantum Poincare recurrences 
with a = 1 scales like S(u) oc 1/uj. A similar behavior of noise has been observed in 



various scientific phenomena ||, for example in the luminosity of stars, the velocity of 
undersea currents, the flow rate of the river Nile, the magnetization of spin glasses and 
the resistance of different solid state devices 0. It is known as 1/ f noise and usually 
extends over several orders of magnitude in frequency, indicating a broad distribution 
of time scales in the system. In the case of quantum Poincare recurrences this property 
stems from the exponentially broad distribution of escape times from some regions of 
the phase space, due to tunneling and localization effects. 

In summary, we have shown that the survival probability for Rydberg atoms in a 
microwave field decays, up to the Heisenberg time tn, in a way similar to the classical 
probability. For t > tu the quantum probability starts to decay slower than the clas- 
sical one, with the exponent of the algebraic decay a — 1. We have given parameter 
values which should allow one to observe quantum Poincare recurrences in microwave 
experiments with Rydberg atoms. Indeed, the thermal beams used with alkali Ryd- 
berg atoms allow one to vary the interaction time by orders of magnitude up to 10 5 
microwave periods pi. Therefore the results of the present paper show that theoreti- 
cal and experimental studies of chaotic Rydberg atoms still represent a challenge for 
fundamental research of quantum chaos. 
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